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Abstract 

A  new  laser  fusion  target  concept  is  presented  with  a  predicted  energy  gain  of  125  using  a  1.3  MJ  KrF  laser. 
This  energy  gain  is  sufficiently  high  for  an  economically  attractive  fusion  reactor.  X-rays  from  high-  and  low-Z 
materials  are  used  in  combination  with  a  low-opacity  ablator  to  spatially  tune  the  isentrope,  thereby  providing 
both  high  fuel  compression  and  a  reduction  of  the  ablative  Rayleigh-Taylor  instability. 


PACS  numbers:  52.58.Ns,  52.25. Nr,  52.35.Py,  52.75. -d 


A  laser  fusion  power  plant  will  probably  require 
energy  gains  from  the  target  of  more  than  100  to 
overcome  the  low  efficiency  of  the  laser.  For 
example,  the  Sombrero  reactor  study'  assumed  an 
energy  gain  of  114,  using  a  7.5%  efficient  3.4  MJ 
krypton-fluoride  (KrF)  laser.  Unfortunately  it  has 
been  difficult  to  design  a  laser  fusion  target  with  these 
high  energy  gains  that  is  at  the  same  time  consistent 
with  the  physics  constraints  of  controlling 
hydrodynamic  instabilities.  We  have  now  developed  a 
target  design  with  a  predicted  energy  gain  of  125 
using  a  1.3  MJ  KrF  laser,  and  preliminary  analysis 
indicates  that  it  is  consistent  with  the  currently- 
understood  constraints  from  the  ablative  Rayleigh- 
Taylor  instability.  Gains  would  be  even  higher  with 
more  laser  energy. 

The  highest  possible  coupling  efficiency  from  a 
laser  to  a  spherically  imploding  DT  fuel  utilizes  the 
“direct-drive”  concept,  in  which  multiple  laser  beams 
directly  and  symmetrically  illuminate  the  spherical 
target.  As  an  example,  the  predicted  laser  absorption 
efficiency  for  our  target  design  is  89%,  and  the 
predicted  rocket  implosion  efficiency  is  10%.  The  net 
coupling  efficiency  of  8.9%  (from  incoming  laser 
energy  to  imploding  DT  kinetic  energy)  compares 
very  favorably  with  the  “indirect-drive”  ignition 
target  design^  for  the  National  Ignition  Facility  which 
has  a  predicted  net  coupling  of  less  than  1.5%. 

In  1995  scientists  at  the  Naval  Research 
Laboratory  (NRL)  completed  a  few-kilojoule  KrF  gas 
laser  called  Nike,  and  demonstrated  that  it  could 
produce  ultra-uniform  laser  illumination^  using  an 
optical  smoothing  technique  that  induces  spatial  and 
temporal  incoherence  in  the  laser  beams.  Each  laser 
beam,  at  1  THz  bandwidth,  has  a  measured  time- 
averaged  intensity  nonuniformity  at  the  focal  plane  of 
only  1%  rms  (excluding  tilt  and  curvature).  With 
multiple  beam  overlap,  the  laser  nonuniformity  drops 


to  approximately  0.2%  rms  (excluding  the  higher¬ 
mode  interference  between  the  laser  beams).  In  the 
succeeding  years'*  Nike  has  been  used  to  ablatively 
accelerate  plastic  foils  on  a  low  isentrope  by  using  a 
pulse  with  a  low  intensity  foot.  With  Nike’s  high 
quality  illumination,  the  mass  nonuniformity 
generated  by  the  laser  imprint  during  the  foot  of  the 
pulse  was  observed  to  be  very  small,  equivalent  to 
less  than  100  A  rms  surface  nonuniformity.  In  a 
spherical  high-gain  implosion  the  imprinted  mass 
perturbations  would  differ  because  of  the  different 
target  materials,  laser  geometry,  and  laser  pulse 
duration,  but  Nike’s  low  mass  imprinting  is 
nonetheless  very  encouraging. 

Since  the  compressed  shell  is  being  accelerated 
inward  by  a  lower-density  plasma  corona,  the  ablation 
interface  between  the  two  regions  is  unstable  to  the 
Rayleigh-Taylor  (RT)  fluid  instability.  This  instability 
can  be  ameliorated^  by  preheating  the  target,  thereby 
reducing  its  density.  The  reduced  density  leads  to  an 
increase  in  the  mass  ablation  velocity,  and  this 
velocity  flow  then  convects  the  RT  instability  away 
from  the  ablation  surface.  A  previous  target  design® 
proposed  using  an  all-DT  target  and  preheating  it  with 
an  enhanced  shock  wave,  but  the  shock  raised  the 
isentrope  of  both  the  ablator  and  the  fuel.  The 
increase  in  the  fuel  isentrope  then  led  to  a  reduced 
spark-plug  compression  and  a  reduced  energy  gain. 
There  appeared  to  be  a  basic  conflict  between  the 
requirement  to  preheat  the  ablator  and  the 
requirement  to  not  preheat  the  main  fuel  region. 

Recently  we  suggested  [Ref.  6]  that  one  way  of 
resolving  this  conflict  would  be  to  spatially  tune  the 
target  isentrope  with  x-rays  instead  of  shocks,  thereby 
putting  part  of  the  ablator  on  a  higher  isentrope  than 
the  main  DT  fuel.  This  proposal  has  now  been 
developed  into  the  design  and  calculations  of  a 
specific  target  to  be  presented  here,  and  illustrated  in 
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Fig.  la.  The  DT  fuel  is  surrounded  by  an  ablator  that 
is  a  CH  foam  (~  10  mg/cm^)  filled  with  frozen  DT; 
the  ablator  can  be  written  chemically  as  CH(DT)g4. 
The  ablator  is  surrounded  by  a  thin  plastic  shell  to 
confine  the  DT,  and  then  by  a  thinner  high-Z 
overcoating  such  as  gold. 

The  CH(DT)g4  ablator,  mostly  hydrogen  isotopes, 
provides  the  low  opacity  that  is  needed  for 
penetration  by  soft  x-rays.  During  the  first  few 
nanoseconds  of  the  laser  pulse,  the  gold  overcoat 
heats  to  ~  70  eV,  producing  broadband  x-radiation 
that  penetrates  into  the  ablator  below  the  K-edge  of 
cold  carbon  which  is  at  285  eV.  The  gold  is  blown  far 
from  the  target  during  the  foot  of  the  pulse  and  does 
not  produce  significant  x-radiation  after  ~  10  ns.  We 
found  that  300-350  A  of  gold  was  needed  for 
significant  ablator  preheating  and  RT  reduction;  a  200 
A  overcoat  produced  negligible  ablator  preheat;  a  400 
A  overcoating  produced  enough  x-rays  to  heat 
through  to  the  fuel  layer  and  prevent  ignition.  In 
addition  to  x-ray  preheating,  the  gold  increases  the 
early-time  separation  between  the  absorption  and 
ablation  surfaces  from  <  10  |im  to  ~  200  |im.  This 
separation  should  relax  the  requirements  on  the  laser 
uniformity.  The  gold  has  the  side  benefit  of  increasing 


300  AAu@  19.3  g/cc 
1  am  CH@  1.07  g/cc 

CH(DT)g4; 

@0.25775  g/cc 
DT  ice  @  0.25  g/cc 

DT  gas  @  0.3  mg/cc 

(a) 


195 1  am 
1690  am 
1500am 


Fig.  1  (a)  High  gain  direct-drive  target  parameters  using  a  1 .3  MJ 

KrF  laser.  The  ablator  is  a  CH  foam  at -10  mg/cm^  filled  with 
frozen  DT.  (b)  Laser  pulse  shape.  The  dashed  lines  at  21.2  ns  and 
25  ns  indicate  when  the  laser  beam  is  zoomed  inward;  i.e.,  it  is 
reduced  in  size  to  match  the  imploding  target. 
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the  infrared  albedo  of  the  target  during  its  injection 
into  the  reactor  chamber. 

The  laser  pulse  shape  is  shown  in  Fig.  lb.  To 
enhance  the  laser-target  coupling  efficiency  we 
utilized  an  important  concepf  called  “zooming”:  the 
laser’s  focal  diameter  was  reduced  in  size  to  follow 
the  target  inward,  first  to  74%  and  then  to  54%  of  the 
initial  spot  diameter.  This  zooming  is  conceptually 
easy  to  achieve  with  the  beam  smoothing  technique 
used  on  a  KrF  laser  such  as  Nike,  because  the 
intensity  distribution  at  the  target  is  an  image  of  an 
aperture  located  at  the  low  energy  stage  in  the  laser. 
Zooming  is  achieved  by  using  multiple  apertures, 
with  the  image  of  the  largest  aperture  relayed  first, 
followed  by  successfully  smaller  apertures.  Electro¬ 
optic  techniques  may  eventually  provide  continuous 
zooming.  Zooming  reduces  the  laser  energy 
requirement  for  this  target  from  2.1  MJ  to  1.3  MJ,  and 
raises  the  predicted  gain  from  72  to  125.  (Shock 
breakout  at  the  inner  surface,  and  the  inward 
acceleration,  begin  at  20.5  ns.  The  shell  implosion 
velocity  reaches  3.4xl0^cm/s.  The  average 
acceleration  inward  during  the  main  laser  pulse  is 
5.74x10^®  cm/s^.  Bang  time  is  28.17  ns.  The  hot  spot 
ignition  radius  is  60  |im  at  an  average  density  of  62 
g/cm^  The  maximum  fuel  density  reaches  560  g/cm^, 
and  the  total  pR  of  the  target  is  1.8  g/cm  .  The 
predicted  yield  is  160  MJ.) 

FASTID  (the  spherical  one-dimensional 
computer  model  used  in  the  above  study)  solves  the 
mass,  momentum,  and  energy  fluid  equations  using 
the  FCT  (Flux  Corrected  Transport)  formalism,  with  a 
sliding-zone  Eulerian  mesh  that  provides  good 
resolution  in  the  physically  important  regions^  Easer 
transport  and  absorption  use  ray  tracing  and  inverse 
bremsstrahlung.  Multigroup  diffusion  models  are 
used  for  both  radiation  and  alpha  particle  transport. 
Electron  and  ion  heat  conduction  use  the  Spitzer  form 
with  flux  limiting.  Tabular  lookups  are  used  for  real 
material  equation-of-state  and  radiation  opacity.  The 
STA  model  is  used  to  generate  tabular  ETE  Planck 
and  Rosseland  opacities,  and  then  the  Busquet 
algorithm  is  used®  to  modify  the  data  base  for  non- 
ETE  radiation  physics.  Typically  we  use  45  radiation 
groups,  non-evenly  distributed;  the  groups  are  chosen 
to  resolve  the  dominant  emission  and  absorption 
energy  lines  and  the  continuum. 

Eigure  2  shows  contours  of  constant  isentrope 
during  the  implosion.  The  isentrope  is  approximated 
here  by  a  pressure  ratio  defined  as  a  =  (T. +P,)/Ppd  , 
where  T.  and  P,  are  the  electron  and  ion  pressure,  and 
FfD  is  the  Eermi-degenerate  electron  pressure  of  fully 
ionized  material.  (It  is  assumed  that  the  DT  and  the 
CH  foam  have  been  uniformly  mixed.)  Eig.  2  shows 
that  the  CH(DT)g4  ablator  is  preheated  by  the  soft 
gold  radiation  during  the  first  10  ns  of  the  laser  pulse, 
and  then  by  the  carbon  radiation  during  the  high 
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Fig.  2  Contours  of  constant  isentrope  show  how  the  gold  and  the 
carbon  radiatively  preheat  the  ablator  material.  The  main  fuel 
stays  near  the  minimum  Fermi-degenerate  isentrope.  The  ordi  - 
nate  is  the  integrated  shell  mass  from  the  center  of  the  target. 

intensity  portion  of  the  laser  pulse.  The  pure  DT  fuel 
stays  on  a  very  low  isentrope,  with  a  ~  1.0. 

We  had  originally  designed  the  target  to  rely 
solely  upon  gold  x-rays  to  heat  the  ablator.  However 
parameter  studies  in  which  we  varied  the  carbon 
density  and  gold  thickness  show  that  there  is  also 
significant  heating  of  the  target  by  carbon  x-rays,  and 
the  optimal  preheating  uses  a  combination  of  these 
two  sources.  There  is  about  a  factor  of  two  flexibility 
in  the  choice  of  the  CH  foam  density:  above  ~  20 
mg/cm^  there  is  significantly  less  x-ray  preheating; 
below  ~  5  mg/cm^  the  x-rays  penetrate  to  the  inner 
edge  of  the  ablator  which  then  thermally  preheats  the 
main  DT  fuel.  In  the  base  case,  with  10  mg/cm^,  the 
Attwood  number  at  the  ablator/fuel  interface  is  less 
than  0.1  during  most  of  the  inward  acceleration.  Thus 
the  RT  at  the  ablator/fuel  interface  (a  separate 
instability)  is  probably  not  important  for  this  target. 

The  laser  fusion  community  has  developed  a 
variety  of  methods  for  assessing  the  ablative  RT 
instability.  The  simplest  guideline  is  the  distance- 
moved-over-thickness,  DMOT  =  (Rg  -  R)/AR,  where 
Rg  is  the  initial  outer  radius  of  the  target,  R( t)  is  the 
mass-averaged  shell  radius,  and  AR(t)  is  shell 
thickness  measured  to  the  1/e  values  from  the  peak 
density.  DMOT  is  proportional  to  J  .^kg  dt,  the  number 
of  classical  e-folds  of  the  RT  instability,  for  a 
perturbation  wavelength  l/k  that  is  comparable  to  the 
shell  thickness.  Another  commonly-used  parameter  is 
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the  in-flight-aspect-ratio,  IFAR  =  R/AR  which  is 
closely  related  to  J  ^kg  dt .  These  two  parameters  are 
plotted  in  Fig.  3.  At  24  ns,  near  the  time  of  peak 
inward  acceleration,  the  DMOT  reaches  40  and  the 
peak  IFAR  reaches  60.  These  values  then  decrease  as 
AR  thickens  due  to  radial  convergence.  Whether  the 
values  are  acceptable  for  direct-drive  targets  depends 
in  practice  upon  the  level  of  nonuniformities  in  the 
laser  and  in  the  target  fabrication. 

The  next  level  of  accuracy  and  complexity  uses  a 
dispersion  relation  for  the  RT  growth.  We  used  the 
formula,  y  =  ^kg  /  (l  +  kl)  -  3k[dm/dtj  /  where  g  is 

the  inward  acceleration,  L  is  the  shortest  density  scale 
length  at  the  ablation  front,  k  =  1/ R{t)  where  e  is  the 
mode  number,  and  [dm/dt)  /  is  the  mass  ablation 
velocity  at  the  peak  density.  Weber  et  al'°  found  that 
the  above  formula  is  a  good  fit  to  their  hydrodynamic 
computer  simulations.  Since  the  main  topic  of  this 
article  is  the  ablative  RT  instability,  we  specified  the 
mass  inhomogeneity  at  20  ns,  the  beginning  of  the 
inward  acceleration  phase.  For  a  mode  spectrum, 
including  both  laser  imprinting  and  fabrication,  we 
assumed  that  the  rms  amplitude  of  the  perturbation  in 
Angstroms  is  given  by: 

A(t)=jE(2l+l}(A,r^^^f 

where  i=2",n  =  l,...,ll-  Following  the  suggestion  of 
Haan",  all  modes  within  a  factor  of  two  have  been 
combined  coherently  and  considered  as  a  single 
mode.  Modes  farther  apart  are  assumed  to  have 
independent  phases.  The  A^^  are  chosen  to  be  equal  at 
20  ns,  with  their  weighted  rms  sum  A(t=20)  =  500  A. 
Each  of  the  A^  then  grow  according  to  the  above  RT 
dispersion  relation.  Figure  4  shows  A(t)/AR(t),  for  the 
case  of  A(t=20)  =  500  A.  The  ratio  only  reaches  18% 
at  26.5  ns,  near  the  peak  inward  acceleration,  with  5.8 
e-folds  of  growth.  We  conclude  that  our  target  may 
survive  the  RT  instability,  if  the  total  perturbation 


Fig.  3  In-flight-aspect-ratio  (IFAR)  and  distance-moved-over- 
thickness  (DMOT)  during  the  shell’s  acceleration  inward.  They 
reach  a  maximum  near  the  peak  acceleration  inward. 


Fig.  4  The  ratio  of  the  rms  amplitude  of  the  shell  perturhation  to 
the  shell  thickness,  versus  time,  using  a  model  dispersion  relation 
for  the  Rayleigh-Taylor  instability.  The  amplitude  is  predicted  to 
remain  less  than  the  shell  thickness  if  the  initial  perturbation  is 
less  than  500  A. 

amplitude  is  less  than  500  A  at  the  end  of  the  foot  of 
the  laser  pulse,  and  if  the  initial  mode  spectrum  and 
RT  growth  rate  formulas  are  in  the  final  analysis 
similar  to  what  we  have  assumed.  (The  results  in  Fig. 
4  also  contain  a  model  [11]  to  evaluate  the  effects  of 
saturation,  but  the  predicted  perturbation  amplitudes 
are  low  enough  so  that  the  saturation  effects  are  not 
important.) 

We  also  considered  two  other  dispersion  relations 
for  the  growth  rate  y:  the  curve-fit  formula  of 
Takabe'^  in  which  L  ^  0  and  there  is  a  factor  of  0.9  in 
front  of  the  square  root;  and  the  analytically  derived 
and  more  complex  formula  of  Sanz^^.  The  Sanz 
formula  gave  the  largest  A(t)/AR(t)  (38%)  and  the 
Takabe  relation  the  lowest  (10%).  There  are  of  course 
basic  physics  limitations  in  all  of  these  dispersion 
relations,  so  we  can  only  conclude  that  there  is  some 
reasonable  possibility  that  our  target  design  concept 
will  survive  the  ablative  RT  instability,  if  the  imprint 
in  the  foot  of  the  pulse  can  be  held  to  ~  500A. 

The  next  level  of  accuracy  and  complexity  would 
be  a  2D  (two-dimensional)  arbitrary-Lagrangian- 
Eulerian  code,  or  an  adaptive-grid-Eulerian  code. 
These  codes  can  calculate  the  laser  imprinting  in  the 
foot  of  the  laser  pulse,  follow  the  RT  instability  into 
the  late-time  nonlinear  regime,  and  include  the  multi- 
mode  incoherent  laser  illumination.  We  have  such  an 
adaptive-grid  Eulerian  code  at  NRE  called  EAST2D 
[8].  An  accurate  resolution  of  a  wide  spectrum  of 
modes  is  demanding  in  computer  resources, 
especially  when  it  includes  a  wide  range  of  physics 
packages  and  accurate  radiation  transport.  However 
the  most  difficult  challenge  has  not  been  the  late-time 
nonlinear  evolution  of  the  RT  instability,  but  rather 
the  simultaneous  calculation  of  the  early-time  blowoff 
of  a  300  A  thick  gold  layer  and  the  ablative  RT 
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perturbations  that  has  important  wavelengths  of 
hundreds  of  microns.  NRE  scientists  hope  to  present 
some  2D  computer  analysis  of  this  target  design  in  the 
near  future  using  the  eAST2D  code,  but  there  are  no 
results  yet. 

Using  the  short  laser  wavelength  of  1/4  pm 
provides  the  maximum  absorption  efficiency  and  the 
maximum  rocket  implosion  efficiency,  along  with  the 
lowest  risk  of  laser-plasma  instabilities.  A  KrE  laser 
at  1/4  pm  also  has  the  advantage  of  flexible  zooming 
that  enhances  the  coupling  efficiency  and  gain. 
However  there  is  some  interest  in  the  fusion 
community  in  the  long-term  potential  of  a  1/3  micron 
solid  state  laser.  We  have  therefore  taken  the  same 
target  described  in  Eig.  1(a),  and  retuned  it  for  1/3 
micron  light,  assuming  the  same  two-step  zooming. 
Zooming  of  individual  laser  beams  is  not  practical  for 
solid  state  lasers  because  they  can  not  relay  an  image 
from  the  laser’s  front  end,  but  in  principle  one  could 
triple  the  number  of  laser  beam  lines,  and  use  each 
third  to  produce  different  image  sizes.  The  yield  of 
the  target  with  1/3  pm  and  zooming  was  155  MJ, 
nearly  the  same  as  with  1/4  pm.  However  the 
absorption  and  rocket  efficiencies  were  sharply 
reduced,  and  the  laser  energy  required  to  drive  this 
target  increased  to  1.8  MJ.  The  energy  gain  thus 
dropped  to  78.  To  compensate  for  this  lower  gain  one 
would  need  a  higher  laser  efficiency  and  a  lower 
capital  cost  per  Joule. 

In  this  paper  we  have  put  forward  for  the  first 
time  a  high  gain  target  that  may  also  provide 
sufficient  control  of  the  ablative  RT  hydrodynamic 
instability.  Gains  well  above  100  are  possible  through 
the  use  of  (a)  direct-drive  laser-target  coupling;  (b) 
controlled  levels  of  radiative  preheat  in  a  low-opacity 
ablator,  (c)  a  1/4  pm  laser  wavelength,  and  (d) 
zooming  so  that  the  image  size  follows  the  target  as  it 
implodes.  The  ablative  RT  instability  that  we  have 
addressed  in  this  article  is  important  during  the 
implosion  phase  which  begins  at  20.5  ns,  when  the 
inside  edge  of  the  shell  begins  to  accelerate  inward. 
We  are  still  evaluating  the  first  20.5  ns  of  the  pulse, 
the  so-called  “imprint”  phase,  in  which  the  laser 
nonuniformities  induce  Richtmyer-Meshkov  modes 
and  seed  the  ablative  RT  instability.  As  we  noted 
above,  the  2D  computational  analysis  of  this  imprint 
phase  is  challenging  because  of  widely  differing  scale 
lengths  (a  thin  gold  coating  and  long  laser 
perturbation  wavelengths).  Our  2D  codes  also  assume 
that  the  laser  and  target  perturbations  are  in  the  form 
of  stripes;  a  3D  simulation  will  eventually  be  needed 
to  model  the  real  laser  speckle  and  RT  instability. 

We  have  not  evaluated  here  the  separate  RT 
instability  that  occurs  as  the  dense  shell  decelerates 
against  the  hot  central  ignitor.  Because  of  density 
gradient  stabilization,  this  deceleration  instability  is 
most  important  for  mode  numbers  <10,  and  is 
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largely  decoupled  from  the  ablative  instability  which 
is  most  dangerous  for  mode  numbers  ~  100.  Our  2D 
implosion  studies  suggest  that  the  deceleration  RT  is 
unimportant  if  one  uses  KrF  lasers  since  they  also 
have  good  beam  uniformity  in  the  low  mode 
perturbations.  The  major  uncertainty  with  the 
deceleration  RT  is  whether  one  can  obtain  sufficiently 
accurate  laser  beam  alignment  and  target  fabrication 
sphericity.  Our  computations  show  another  fluid 
instability'^  called  the  “radiating  plasma  structure”, 
which  can  introduce  deleterious  density  waves  and 
preheating.  Some  evidence  of  this  structure  is  visible 
in  Fig.  2  as  the  notch  in  the  a=1.5  isentrope  at  14  ns. 
However  by  carefully  tuning  the  laser  pulse  we  were 
able  to  avoid  the  deleterious  effects  of  this  instability. 

We  also  have  not  discussed  the  various  laser- 
plasma  instabilities.  Preliminary  analysis  indicates 
that  these  instabilities  are  near  threshold  and  may  be 
tolerable  because  of  the  KrF  laser’s  short  1/4  |im 
wavelength  and  the  modest  peak  laser  intensity  (<10'^ 
W/cm^),  but  more  theory  and  experiment  are  required. 
We  have  also  assumed  that  the  CH  foam  and  the  DT 
in  the  ablator  are  uniformly  mixed  from  the  beginning 
of  the  acceleration  phase.  We  have  some  preliminary 
experimental  and  computational  evidence  that  this  is 
not  correct.  We  are  assuming  that  the  laser  pulse 
shape  can  be  retuned  to  correct  for  the  actual  shock 
velocities. 
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